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ABSTRACT 



I prove the recently conjectured relation between the 2 x 2-niatrix differential operator 
L = d^ — U, and a certain non-linear and non-local Poisson bracket algebra (V^-algebra), 
^ ' containing a Virasoro subalgebra, which appeared in the study of a non-abelian Toda field 

theory. In particular, I show that this V^-algebra is precisely given by the second Gelfand- 

• 1—1 . 

/\f ■ Dikii bracket associated with L. The Miura transformation is given which relates the second 

b : 

. 5t , to the first Gelfand-Dikii bracket. The two Gelfand-Dikii brackets are also obtained from 

the associated (integro-) differential equation satisfied by fermion bilinears. The asymptotic 
expansion of the resolvent of (L — ^)\E' = is studied and its coefficients Ri yield an infinite 
sequence of hamiltonians with mutually vanishing Poisson brackets. I recall how this leads to a 
matrix KdV hierarchy which here are flow equations for the three component fields T, V'^, V^ 
of U. For V^ = they reduce to the ordinary KdV hierarchy. The corresponding matrix mKdV 
equations are also given, as well as the relation to the pseudo-differential operator approach. 
Most of the results continue to hold if ?7 is a hermitian n x n-matrix. Conjectures are made 
about n x n-matrix, ?7i* -order differential operators L and associated Vj'^ ,„)-algebras. 
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1. Introduction 

During the last ten years there has been made tremendous progress in understanding the 
close relations between integrable hierarchies of KdV type (see e.g. ref. [1]) and field theory, 
in particular conformal field theory. Early examples include the correspondence between the 
second hamiltonian structure of the KdV equation and the Virasoro algebra [2,3], with the free 
field realization of the latter providing the Miura transformation for the former. (Analogous 
results hold for the supersymmetric extensions, see e.g. [4]). This pattern was generalized with 
the discovery of the relation between Toda field theories, W^-algebras and the m Gelfand- 
Dikii hierarchies (m reduction of KP) [5,6,7]. The general association of a hamiltonian 
structure (Poisson bracket) with a differential operator is mainly due to Gelfand and Dikii 
[8,9] and Drinfel'd and Sokolov [10]. The aim of the present paper is to elaborate on yet 
another generalization of these patterns. 

In the recent study of the non-abelian Toda field theory as a model for strings propagating 
on a black hole background [11], a new Poisson bracket algebra of the conserved currents was 
discovered. It includes the Virasoro algebra with classical central charge as a subalgebra. This 
algebra contains, besides the stress-tensor T, two other spin-two currents V^ and V^ . It is a 
non-linear and non-local algebra and was called l^-algebra to emphasize the similarities (non- 
linearity) and differences (non-locality) with the well-known W^-algebras [12]. The algebra 
reads: 

l-HT{a) , T{a')} = (9. - d^,) [T{a')6{a - a')] - ^-6"'{a - a') 

7-2{T(a) , V^a')} = {d, - d,,) [V^{a')6{a - a')] 

7-2{\/±(a) , V^{a')} = e{a - a')V^{a)V^{a') (1.1) 

7-2{\/±(a) , V^{a')] = -e{a - a')V^{a)V^{a') 

+ {d, - d,>) [na')6{a - a')] - ^6'" {a - a') . 

Here 7^ is a (classically arbitrary) scale factor related to the classical central charge c by 
c = 127r7^ , as can be easily seen by looking at the modes L„ = 7" /^^ do'[T(r, a) + j\e^'^^'^^"'. 
To be definite, the coordinate a is supposed to take values on the unit circle S\ and t{a — a') 
is defined to be the unique periodic, antisymmetric and translationally invariant function 
obeying dada'e{a - a') = 25'{a - a'), namely e{a - a') = 9{a - a') - e{a' - a) - 2^ 



= ^^m^n:^^^^ ^ ■ However, for many of the results presented here, a E S^ is not 
essential, and can be replaced by cr G R. 

In ref. [11] a free field realization of the T and V was given (see section 3 below). Using 
this realization it was shown that simple (vertex-operator type) exponentials ifji and ■?/'2 of the 
free fields satisfy a matrix differential equation 



()i 





:i.2) 



Below, I will often use the matrix notation 



L^ = 0, L = d^-U , U 



T 

-V2V- 




^ 




:i.3) 



where, if not indicated otherwise, 8 = 8^, and U and \1' depend on a (as well as on other 
parameters ti,t2, ■ ■ ■ which will not be relevant until section 6). The dependence on r, inherited 
from the conformal field theory is trivial since U and ^ only depend on the combination a + r. 
Actually, as usual (see below), this is the first flow of the matrix KdV hierarchy: -nj-U = daU, 
so that r is identified with ti. 

The plan of this paper is fairly well described by the abstract. 

2. The Gelfand-Dikii brackets associated with L = d^ — U 



In this section, following refs. [9,10,6,13], I define the first and second Gelfand-Dikii brack- 
ets (hamiltonian structures) associated with the 2 x 2-matrix second-order differential operator 
L. Then I show that the second Gelfand-Dikii bracket is identical to the V^-algebra (1.1). 

Let / and g be polynomial functionals on the space of second-order differential operators 
L, i.e. polynomial functionals of U (and its derivatives). One defines the pseudo-differential 



* To be precise, the equation of ref. [11] is a differential equation in u = t + a. However, for the present 
purpose, we can consider the "01) V'2 to be exponentials of the left- moving (i.e. u-dependent) parts of the 
free fields only, so that daipi(u) = duipi(u), and the equation given here is equivalent to that given in ref. 
[11]. 



operator 



Xf = d-'^Xi + 3-^X2 



X -^f 



(2.1) 



where d ^d = dd ^ = 1 and -Aj is defined as 



su 



1 5 


1 5 


2 ST 


V2SV- 


1 S 


1 3 


V2^V+ 


2W 



(2.2) 



t 



so that jjj Jtr f/"" = nU"^^^, and X2 is determined (cf. [6,9,13]) by requiring res[L, Xj] = 0. 
As usual, the residue of a pseudo-differential operator, denoted res, is the coefficient of d^^. 
One has 



x; ^ 1 (U. 



2\5U 



1 

+ 2 



' 6U 



(2.3) 



Integrating this equation yields X2. Here, one observes a new feature as compared to the 



scalar case: since in general 



^' su 



7^ 0, X2 will be given by a non-local expression involving 



an integral. This is the origin of the non-local terms (~ t{a — a')) in the K-algebra (1.1). 
The first Gelfand-Dikii bracket is defined by 



{/,fl'}GDi = a / datr res {[L,Xf]+Xg) 



(2.4) 



where [■, ■] is the commutator and the subscript + indicates, as usual, to take only the differen- 
tial operator part (no negative powers of d). The parameter a is a real, non-zero scale factor, 
to be fixed later on. It is easy to obtain 



{f.9] 



GDI 



-2a I datr gd^ 



(2.5) 



I assume that one can freely integrate by parts, which is true in particular if o" G S*^. If o" G R, 
one has to require that U vanishes sufficiently fast (e.g. exponentially) as o" ^ ±00. 



f On the circle S^ e.g. 9 ^ is well defined on functions without constant Fourier mode, i.e. 



/(^) = Em^o/' 

-E 



m^O ' 






One easily sees that {d ^f)((7) = /dcr' ^^{<^ ^ '^')I{'^') with e{a — a') 



I This condition is necessary since the coefficient of 9 in i vanishes. 
§ using the well-known formula d^^f = E^o(~)"/'^"''^^^ " 



The second Gelfand-Dikii bracket is usually more interesting. In analogy with the standard 
procedure [9,10,6,13] I define it as follows. 



{f,9}GD2 = aJdatT Tes{L{XfL)+Xg-{LXf)+LXg) . 
It is straightforward to obtain 



{/,.}oo.^a/d.tr|lf.3^ + i 



' 6U 



d-' 



'6U 






(2.6) 



(2.7) 



where the (9 ^ is meant to act only on 



U ^ 



Inserting the definitions of the 2 x 2-matrices U and -^ one obtains 

r n ^ f T f ^f r^Sg ^f r. Sg ^f r. Sg 

l/,fl'|GDl = -a I da \ -T^d^:^ + tttt^tt^ + -tt^o- 



6T 5T 6V+ 5V- 5V' 5V+ 



and 



a f 
{f,g}GB2 = ^2 / ^^^^ 



_ 15/^3 5^ _ 15/ ^3 5g 



-d' 



+ T 



2 6T 6T 2 5V+ 6V- 
Sfjg , Sf ^ 6g 



1 Sf_gs Sg 



2 5V- 6V+ 



:d 



d- 



6T 5T 6V+ 6V 



'f-a-'' 



6V- 6V+ 



if^g) 



+ v+i^-ld ^^ ■ ^^ -^^ 



+ v- 



5T 5V+ 5V+ 6T 



d— -U ^ g) 



ST 5V- 







dada'ela - a) V+ 



5V+ 



V- 



5V- 



MW 



r+ 



l9_ 

6V+ 



V- 



-Jg_ 
6V- 



{21 



o- 



(2.9) 



Then, taking /, g to be T, V^ or V one has the 

Proposition 2: The second Gelfand-Dikii bracket (2.9) coincides with the l^-algebra (1.1) 
provided one chooses 



a= -27 



-2 



(2.10) 



Henceforth I will adopt this choice and the only Poisson bracket used is the second Gelfand- 
Dikii bracket, unless otherwise stated. 



3. Free field realization and the Miura transformation 

In ref. [11] the l^-algebra (1.1) was obtained from a field tlieoretic realization. After a 
series of field redefinitions, which need not concern us here, the following free field realization 
of the generators was obtained 

^ z-i v2 (3,1) 

\/± = i(v^v^'3 - d) [e^'^^^(v^; T ^V^'2)" 

where, for the present purpose, d = do-, ^'i = dipi and (fi = (/9i(r + a) is a chiral half of a free 
field. Then it was shown in ref. [11], and can easily be verified again, that the transformation 
(3.1) maps the following free-field Poisson brackets to the brackets (1.1): 

Wiia), ^j{(r')} = -—6^je{a - a') . (3.2) 

Note that this implies 

Wi{(^),Vji(^')} = I'^SijS'ia - a) (3.3) 

or, if a, a' E S^, i{(pn,(p'^} = n5^^6n+m,o for the Fourier moded defined by (p'j{o-) = 
7(27r)~^'^ ^j^y9ne~*"'^. These are harmonic oscillator Poisson brackets. Of course, (3.3) has 
the structure of the first Gelfand-Dikii bracket (2.8). Thus (3.1) maps the first Gelfand-Dikii 
bracket to the second one, or, in other words: 

Proposition 3: The transformation (3.1) provides a Miura transformation for the two hamil- 

tonian structures associated with L = d'^ — U. 

The Miura transformation (3.1) can be nicely written in matrix form. Define 

V2 \ V2p- -p3 I (3 4) 

P3 = <^3 > V± = ^e^'^^^'(<^i T iv'2) ■ 




-k Of course, on could replace d^ by 9„ — \{da + dr) everywhere and consider (pi to be the full free field, 
i.e. a sum of a right- and a left-moving part. This is the approach of ref. [11]. 



Then 



U = P'^ + P' 



(3.5) 



and the Poisson brackets of the p are 



7"^{P3(f^),P±(f^')} = 
l^'^{p±{o-),p±{cr')] = e{a - a')p±{a)p±{a') 

1^^{p±{(^),Pt{(^')} = -e(^ - ^')p±(^)Pt(^') 



(3.6) 



6' (a - a') 



.:..J 



The matrix Miura transformation (3.5) allows us to to rewrite the differential operator L as 



L = 92-f/= (d + P){d-P) . 



(3.7) 



Note that (3.1) implies that V^ and V^ are complex conjugate if the ipj are considered 
real. All equations are compatible with this assumption. In all of the following, I will suppose 
that {y^Y = V^ . In particular, this implies that f/ is a hermitian matrix, so that L is a 
hermitian operator. Let me also note that one solution of the equation L^ = is easily 
obtained. If T and V"^ are given by (3.1), a solution is [11] 



ipi = exp 
■02 = exp 



— =(V9l -?^2-^3) 



[V2 



{ipi +iip2 - V^3) 



(3. 



ri 



t Strictly speaking, as a mapping between the two Gelfand-Dikii brackets, the Miura transformation 
is the transformation (3.1) from the f'^ to the T,V^. Since, due to the e{a — cr'j-terms, the Poisson 
brackets (3.6) of the p are not exactly the first Gelfand-Dikii bracket, equation (3.5) is not really a Miura 
transformation. However, by abuse of language, it is convenient to refer to it nevertheless as the matrix 
Miura transformation. 



4. The integro-differential equation for (fermion) bilinears 



It is often useful to introduce a spectral parameter A and consider the eigenvalue problem 



-L^ = \^ 



(4.1) 



(with A G R, since f/ is a hermitian matrix). It is suggestive to think of this as a Schrodinger 
equation for a fermion, and thus natural to consider also fermion bilinears. Define 



XI = (V^i>l + V^2V'2) = ^+^ 



(4.2) 



The reason to choose the Pauli matrices o"+, c- and the unit matrix (rather than as) lies in 
the form of the matrix U = T — \P1gJV^ — \f2oJV^. Let x be the vector with components 
Xi; X2, X3- Oiie then has the following 



Lemma 4: Assuming that {V 



+\* 



V , X satisfies the following integro-differential equation 



Ax = AiVx , N = 2 



n o\ 

1 
\0 1 0/ 



d 



(4.3) 



where A is the operator 



t 




+ 



I T 

V+ 

\v 




T 



T 
/ 



d 



(4.4) 



\ As usual, d and d ^ are meant to act on everything on their right in Ax. 
{V-d-^V-X2){<y)=V-{<j)id<j'\e{a-a')V-{a')x2{<j'). 



For exemple, 



The matrix elements of this operator are simply the second Gelfand-Dikii brackets: 



((tIAIo- ) = 7 



/ {T{a\T{a')] {T{a\V+{a')] {T{a\V-{a')] \ 

{V+{a),T{a')} {V+{a),V+{a')] {V+{a),V-{a')] 
\{V-{a\T{a')] {V-{a\V+{a')] {V-{a\V-{a')] ) 



(4.5) 



while the matrix elements {(j\N\a') obviously are the first Gelfand-Dikii brackets. 

Proof: Equation (4.5) follows obviously from (4.4) using (1.1). To show that (4.3) is satisfied 
is equally straightforward. For example, one of the equations is 

[^3 _ 29(T - A) - 2(T - \)d]i)li)2 + V2{dV- + \/-9)(V^i>i + V^^V^a) 

(4.6) 

= -AV-d-^{v^i)li)2 - v-i)*2iJi) 

Using (4.1) in component form, i.e. [d'^-{T-X)\i)i+^/2V+i)2 = 0, [d'^-{T-X)]ilj2+V2V-^i = 
0, and the complex conjugate equations (with (V^)* = V^) one easily sees that the l.h.s. of 
(4.6) equals ^/2V^ {^^d^i —d^i^pi — ijj2'^^2 + d^2^2) ■ This equals the r.h.s. of (4.6) since one 
has \/2d{tljldtlJi — d^pltpi — il'2^'4'2 + <9'?/'2'^2) = AV'il)2'4'i — W^il)\il)2- The two other equation 
are proven in the same way. 

Let me emphasize again that the integro-differential equation for the fermion bilinears has 
reproduced the first (A^) and second (A) Gelfand-Dikii brackets. This is not too surprising in 
view of the important role played by fermion bilinears in the usual KP hierarchy [14]. 



5. Asymptotic expansion of the resolvent 

This section is inspired by the classical work of Gelfand and Dikii [8] where the asymptotic 
expansion of the resolvent is given for scalar U'. Most results of [8] have a straightforward 
generalisation to the present matrix case, although some of the proofs have to be modified. The 
matrix case has been studied by Olmedilla, Martinez Alonso and Guil [15] who, in particular, 
were the first to prove Proposition 5.2 below. Since our proof is very compact, and in order 



I use "scalar" as opposed to "matrix" . Of course, under a coordinate transformation, it is most convenient 
to view U as transforming as a quadratic differential with an additional Schwarzian derivative term. 



to keep this paper self-contained, it will be presented here. For a hermitian n x n-matrix U 
define the resolvent R as 



R{x,y-i) = {x\{-d' + U + ir'\y) (5.1) 



which is a solution of 



{-dl + U{x) + i)R{x, y- = 6{x - y) . (5.2) 



R{x^ y; ^) is well-defined for any ^ such that — ^ is not an eigenvalue \a of —d + U. If the 
spectrum is discrete, as I will suppose, R{x, y; ^) is a meromorphic function in the complex ^- 
plane. Let Rq be the resolvent for f/ = 0, e.g. for x,y eH one has Roix, y; ^) = y7?^^ ^ 
Then R has the (formal) expansion R = RoYlT=oi~) i^^o) ■ "^^^ k = 1 term e.g. is 
— J dzRQ{x, z; ^)U{z)Ro{z,y; ^), and developping U{z) = ^ ^"^-j^U^ ' (x) one finds, for x = y, 
the sum — Yl'n=l Antn+i/V ■ ^^^ can proceed in the same way with the terms k > 1. For other 
boundary conditions, e.g. x,y G S^, similar considerations apply. As a result one obtains the 

Lemma 5.1: The restriction of the resolvent to the diagonal, R{x;^) = R{x,x;C,) has an 
asymptotic expansion for ,^ ^ oo of the form 



CX) 



Rn[u\ 



R^^-^0 = 1^t£^,- (5.3) 

n=0 ^ 

This equation is to be understood as an equality of the asymptotic expansions in powers of 
,^~^/^, disregarding any terms that vanish exponentially fast as ,^ ^ oo. (These exponentially 
vanishing terms do depend on the specific boundary conditions, while the coefficients of the 
asymptotic expansion do not [8].) The coefficients Rn[u\ = Rn{x) are (izjferentm/ polynomials 
in U (i.e. polynomials in U and its derivatives). It also follows from the above construction 
that terms like e.g. U'^U" appear in a symmetrized form ■^(U'^U" + U"U'^). In particular, if U 
is a hermitian 2 x 2-matrix with Un = U22, the same is true for R{x; ^) and Rn{x). 



% For the case x, y G [— tt, tt] with periodic boundary conditions one has Ro(x,y;S_) ~ —^ — 7= 

and one finds for the fc = 1 term —J2'^=i Anfn+i/V + 0{e^y^^), which coincides with the result for 
a;, y G R up to terms that vanish exponentiaUy fast as ^ ^ 00. 



It is well-known from standard quantum mechanics that if {^^(a;)} is a complete set of 
orthonormal eigenf unctions of the eigenvalue problem 

{-dl + U{x))'^a{x) = Xa'^aix) , (5.4) 

i.e. satisfying J2a ^0(^)^0: (^) — ^(^ '~ v)^ then the resolvent has the spectral decomposition 

«(,.,;^).^M£MM, (5,5) 

a ^ ^ 

Note that it follows from the spectral decomposition and the normalization of the '^a. that 

jtTR{x;0dx = J2^^- (5-6) 

Proposition 5.2: For an x n hermitian matrix U let the resolvent R{x, y; $,) be a solution of 
(5.2). Then R = R{x; ^) = R{x, x; $,) satisfies 

R'" - 2{UR' + R'U) - {U'R + RU') + [U, d^\U, R]] = 4^R' (5.7) 

(where R' = dxR{x; C,) etc.) and the coefficients Rn of the asymptotic expansion (5.3) satisfy 

4<+i = C - 2{UR'^ + R'^U) - {U'Rn + RnU') + [f/, d'^U, Rn]] ■ (5.8) 

Proof. It is straightforward to see, using (5.2), that -R(x, y; S,) satisfies the following differential 
equation in x and y: 

{dx + dyfR{x, y; - 2U{x){d, + dy)R{x, y; - 2((9, + dy)R{x, y; 0)U{y) 

- U'{x)R{x, y; - R{x, y; OU\y) + U{x) {d, - dy)R{x, y; (5.9) 

- ((9, - dy)Rix, y- i))U{y) = A^d, + dy)R{x, y; ■ 

Note that all 6' {x — y) -terras cancel. Next, since {dx + dy){dx — dy)R(x,y]^) = U(x)R(x,y;^) — 
R{x,y; OUiy) onehas {dx-dy)R{x,y; ^) = {dx+dyy^{U{x)R{x,y;^)-R{x,y;C)U{y)). Insert- 
ing this into eq. (5.9), taking x = y and observing that (^{dx + dy)^R{x, y; ^)) \x^y = d^R{x; ^) 
one obtains the desired equation (5.7). Equation (5.8) then is an immediate consequence of 
the expansion (5.3). 



* Recall that I suppose a discrete spectrum {Aq}. This avoids, among other things, unnormalizable 
eigenfunctions. However, since most of the results involve only formal algebra in U and its derivatives, 
this is not really essential (cf. ref. [8]). 

10 



Although not obvious from the recurrence relation (5.8), one is garanteed by Lemma 5.1 
that this can be integrated to yield all Rn as differential polynomials in U . The first few are 
(for U a hermitian n x n-matrix) 

R^=-\u 
io 



6U{x) 2 



6U{x] 



JtTR{y;Ody=^Rix;^) 



(5.10) 



i?3 = - — (lOf/3 - 5UU" - 5U"U - 5U'^ + f/(^)) 
64 

i?4 =^(35f/4 _ 21[/2f/" - 21f/"[/ - 28UU"U - 28[/'V - 28UU'^ - UU'UU' 
+ 7f/f/(^) + 7f/(^)f/ + 14f/'f/'" + IAU"'U' + 21f/"^ + f/(6)) . 



Lemma 5.3: Let f/ be a general hermitian n x n-matrix. Then 

A 1 /■ 2n — 1 

- / tr i?n(l/)di/ = ^i?n-l • (5.11) 



Proof: From (5.6) one has ^ / tr R{y- i)dy = ^ E« eTAl = - Ea 5U(x){i+\^y ^tan- 
dard first order perturbation theory gives 5Aq, = J ^J(5f/\1'Q, hence 

^ = ^aK (5.12) 

where {S/SU)ij = 6/6Uji. Now, from (5.5) it then follows that 



(5.13) 



Inserting the expansion (5.3) gives the desired result. Note that when U is not a general 
hermitian matrix, some care has to be taken when defining {6/6U). For exemple, for the 
2 X 2-matrix U of eq. (1.3), {6/6U) is defined by (2.2), and it is easy to see that eq. (5.13) 
continues to hold. 



** Indeed, the last term, related to the non-local terms in the V^- algebra, and caused by the non- 
conimutativity of the matrices, makes it even less obvious than in the scalar case. 

11 



Now let me return to the special case where U and and hence also R{x; ^) both are her- 
mitian 2 x 2-matrices with equal (1, 1) and (2,2) elements. It is then convenient to define a 
correspondence with 3-component vectors as 



M 




b 



(5.14) 



This definition, together with the spectral decomposition (5.5), as well as the property 
-R(2;;0ii = -R(^;022! and the definition (4.2) of x iXi = ^a^a ^^c) imply 

Then we have the following 

Corollary 5.4: It is a simple matter of writing out the components explicitly to see that, if U is 
the familiar 2x2 matrix, eqs. (5.7) and (5.8) are equivalent to the following integro-differential 
equation 

AR{x;O+^NR{x;O = (5.16) 

where A and A^ are defined in eqs. (4.3) and (4.4), and for the coefficients Rn{x) 

ARn + NRn+1 = . (5.17) 

6. The infinite sequence of hamiltonians 

In this section, I will construct an infinite sequence of hamiltonians Hn with {Hn, Hm} = 0, 
see also [15]. (For a discrete, i.e. difference version of the matrix differential operator L an 
infinite sequence of conserved quantities was constructed in [16].) The definitions (2.5) and 
(2.7) of the first and second Gelfand-Dikii brackets are sufficiently general to allow for arbitrary 
hermitian n x n-matrices U. It is also obvious how to generalise the Proposition 6.3 below to 
this case. For simplicity of presentation, however, here I will restrict myself again to the case 
where U is the hermitian 2 x 2-matrix of eq. (1.3). 

12 



That such a sequence of hamiltonians with {Hn, Hm} = exists follows from the next 

Lemma 6.1: The second Gelfand-Dikii bracket of any two spectral parameters A (eigenvalues 
of —d + U) vanishes: 

{Xa,X/3} = 0. (6.1) 

Proof. If Xa = Xp the bracket vanishes trivially. So let Aq 7^ A^. One has by (5.12) and the 
definition (4.5) of A 



r\KAp] = I dxdy {x''{x)f{x\A\y)x^{y) 

(6.2) 
dx {x^{x)f{Axf'){x) . 



By lemma 4, eq. (4.3), this equals A^ J dx {x"{x))'^{^X ){x)- It follows from the antisymme- 
try of the Poisson brackets (or directly from (4.4)) that A is an antisymmetric operator, as is 
obviously also A^. Thus 7~^{Aa, A^} also equals Xa J dx {x°'{x))'^{Nx^){x)- For Aq 7^ A^ this 
implies j dx (x"(a;))-^(A^x^)(a;) = and hence {Aq,, A^} = 0. 

Corollary 6.2: The first Gelfand-Dikii bracket of any two spectral parameters A vanishes, 
too. 

Proof-. This bracket is given by J da; {x^{x)Y'{Nx ){x) which was just shown to vanish. 

Thus the infinite set of Aq is in involution. However, their dependence on U is not very 
useful for practical purposes. On the other hand, the operator-trace of the resolvent operator, 
J da; tr R{x] ^) = j dx tr (a;|(— 9^ + U + ^)~^\x)^ provides a very convenient functional of the 
eigenvalues Aq, only (cf. eq. (5.6)). It follows from the previous lemma that 

{j dx tr R{x- Ojdy tr R{y; ^')} = (6.3) 



and if one defines 



Hn = ^^y^ / dx tr Rnix) (6.4) 



one has, upon inserting the expansion (5.3) into (6.3), the 
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Proposition 6.3: The Hn are in involution: 

{Hn,H^} = Q. (6.5) 

Since tliis proposition is important, I will prove it once again in a different way, along the 
lines of [15]. By Lemma 5.3 of the previous section, eq. (5.11), one has 

w^ = (-^)"«"(^> ■ («■«> 

From Corollary 5.4 one then gets 

Lemma 6.4: The Hn satisfy the recursion relation 

or in components 

2{T,Hn} = d^-^ , 2{V-,Hn} = d^-^ , 2{V\Hn} = d^-^. (6.8) 

Thus, since A defines the second Gelfand-Dikii bracket, and since A and A^ are antisymmetric 
operators, one has 

6Hrt ^ 6Hm 1 f SHn 6Hm+l 



' ^ ' ^ J 5U 5U A J 



SU SU 



4/^^wP^-/^^nr (-' 



f 5Hn-l . SHm+1 -2 r rj u ^ 



Iterating this m times one arrives at 7 {Hq, Hn+m} which vanishes since Hq is constant 
(independent of U). Hence {Hn,Hm} = is proven again. The equalities of eq. (6.9) also 
show that {Hn,Hm} = {Hn,Hm}GD2 = {{Hn, Hm+ijcDi i-G- the Hn are in involution with 
respect to both Gelfand-Dikii brackets. More generally, the above Lemma 6.4 shows that for 
any functional f[U] one has 

4{/, Hn}GD2 = {/, Hn+l}Gm ■ (6.10) 

The model is bihamiltonian. 
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The first few Hn (most easily obtained from (6.4) and (5.10)) are 
H, = UtrU=fT 



H2 = ^ fti U^ = [{T^ + 2V+V-) 

H3 = - I tT (2f/^ + f/'2) = /(2T^ + 12TV+V~ + T'^ + 2V+'V' 

H^ = - I ii (5f/^ + 10f/f/'2 + f/"2) 

= M5T^ + 20y+ V"^ + 60TV+y" + lOTT'' 



(6.11) 



n/2 

+ 20T\/+V"' + 20TV+\/"' + 20TV+V" + T"^ + 2\/+'V"") . 

Remark 6.5: The proof of {Hn,Hm} = is based on the recursion relation (6.7) inherited 
from that of the coefficients Rn of the resolvent. The crucial point is that A and A^ give 
the second and first Gelfand-Dikii bracket. Comparing the recursion (5.8) for the R^ for 
general (hermitian) U with the first and second Gelfand-Dikii brackets (2.5) and (2.7) for 
general (hermitian) U it is obvious that the proof of {H^, Hm} = can be straightforwardly 
generalised to f/ a hermitian n x n-matrix. 

Remark 6.6: In the scalar case (f/ a 1 x 1-matrix) one can show [3] that, if px is a solu- 
tion of the Riccati equation U — \ = p\ + p'^, then the quantity A;(A) = exp(Jj9;^(a;)da;) is 
an eigenvalue of the monodromy matrix of the wave function, and thus only depends on A. 
Hence by the above Lemma 6.1, {log A;(A), log A;(A')} = 0. From the Riccati equation one has 
PX = V^ + E^=i (^^ with c.(-+i) = -(c.(H)' - j:^J^'Jr)^ir-r)^Ji) = u. Hence 
{/cj(""), /cj(")} = 0. It turns out that /cu^^m) ^ g while (-)"+! /cj(2"-i) gives the hamilto- 
nian Hn- I did not succeed to give the same type of proof in the matrix case. However, it is 
easy to verify for n = 1, 2, 3 and 4 the following 

Conjecture 6.7: Let f/ be a hermitian n x n-matrix. Define Vt^^' = U and ^2*'"^+^) = 

_(f]My _ ^m-l ^(r)^(m-r)^ Let Hn = (-)"+4 /tr fi(2n-l)^ rj.^^^ ^^ ^ ^^ 

A natural question is whether the Hn are the only hamiltonians and whether this set is 
complete. Preliminary considerations of some simple other functionals of U, i.e. of T, V^ and 
V^ seem to indicate that this is indeed the case. 
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7. Hierarchies of matrix KdV and matrix mKdV flows 
As usual one defines an infinite hierarchy of flows by 

f^7-W//.}^ (7.1) 

Since {Hr, Hs} = it follows from the Jacobi identity that all flows commute. Also, as usual, 
the flow in ti is trivial: dU/dti = daU, implying that f/ is a function of a + ti (and t2, ts, . . .) 
only. Thus ti can be identified with the "world-sheet time" r. The flow in t2 gives the matrix 
generalisation of the KdV equation: 

1^ = (3f/2 - U'y (7.2) 

or in components 

^ = {3T^-T" + 6V+V-y , ^ = {6TV^-V^"y (7.3) 

OT2 Ot2 

Just as the Virasoro algebra is a subalgebra of the V-algebra (1.1), the KdV equation is simply 
obtained by setting V"^ = 0. Note that since all Hn are symmetric in V'^ and V^ any non-local 
term (~ e{a — a')) that might appear cancels in all flow equations, and the latter are always 
partial differential equations. 

Using the Miura transformation of section 3, eq. (3.5), one also has the hierarchy of matrix 
mKdV flows ^ = 7~^{P, Hr}., in particular 

or in components 

_ = (p3 + 6p3P+P- - p'i) , -Q^ = {3pIp± + 2pip^ - p^) . (7.5) 



-k Matrix KdV equations were already considered a long time ago by Calogero and Degasperis [17]. 

f As before, {•, •} is meant to be the second Gelfand-Dikii bracket. But since A§^ — A^~^{U, Hr\ gt)2 = 

7^^{f7, i/r+ilcDi, one can also use the first Gelfand-Dikii bracket and the next higher hamiltonian 

instead. 
I This also follows from the equivalence with the first Gelfand-Dikii bracket which is local, see previous 

footnote. 

16 



8. Relation with the pseudo-differential operator approach 

As an alternative method to study the hierarchy of matrix KdV flows one can use the 
pseudo-differential operator method [18]. In the scalar case this is very standard technology 
(see e.g. [19]), while here some attention has to be paid to the non-commutativity of U,U' 
etc. I will only sketch how to obtain the matrix KdV equation (flow in t2). Starting with 
L = d'^ — U one deflnes the square root 

l'^ =d- -Ud-^ + -U'd-^ - -(U^ + U")d-^ + . . . (8.1) 

2 4 8 



so that 



and 



33 3^^^ 3^^; 



(L-^)^ = d'--Ud--U' (8.2) 



[(Li)+, L] = -iWU' + WU - U'") = i(3f/2 _ jj^y ^g g^ 



4^ ' 4 

so that the matrix KdV equation reads as expected 

IdL 



ASt, "^^'-^l- <"' 



9. Conclusions and generalisations 

In this paper, I have shown that the recently discovered non-linear and non-local K-algebra 
(1.1) is obtained as the second Gelfand-Dikii hamiltonian structure based on the second-order 
2 X 2-matrix differential operator L = d"^ — U . I observed that the non-locality of the V-algebra 
is a direct consequence of the non-commutativity of matrices. The Miura transformation is 
given by the free fleld representation and relates the flrst and second Gelfand-Dikii brackets. 
It is shown that —L^ = A\E' implies an integro-differential equation for bilinears in ^. This 
latter equation naturally reproduces the two Gelfand-Dikii brackets. Then the asymptotic 
expansion of the (matrix) resolvent is studied, and the integro-differential equation for the 
bilinears translates into a recursion relation for the coefficients Rn- Defining an infinite series 
of hamiltonians as integrals of matrix-traces of the Rn, the recursion relation of the latter 
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becomes a recursion for the hamiltonians. This immediately imphes that the infinite sequence 
of hamiltonians is in involution. The corresponding fiows all commute and give the matrix 
generalisations of the KdV (and via the Miura transformation also of the mKdV) hierarchy. 
Finally, I outlined the connection with the matrix pseudo-differential operator approach of ref. 

[18]. 

There are obvious generalisations of the previous results. The generalisation to the case 
where f/ is a hermitian n x n-matrix was already mentioned at several occasions. More gen- 
erally, one can study the n x n-matrix, -m* -order differential operator L = d"^ — Ufo-xd'^^'^ — 
^(3)'^'^"^ — ... — f/(^), maybe with certain constraints on the form of the matrices Uny (The 
matrix U in this paper was constrained by tr o"3f/ = 0.) The corresponding Gelfand-Dikii alge- 
bra could be called V?^ „j-)-algebra, so that the above V-algebra would be V(2,2)? while the usual 
IVjV-algebras are V(\m- There has been recently very interesting work by Kac and van de Leur 
[18] on the n-component KP-hierarchy and the corresponding m reductions, mainly (but not 
exclusively) in the pseudo-differential operator formalism. The present approach, which is 
more field theoretic and insists on the V-algebraic aspects and the hamiltonian structures, 
should be viewed as complementary. 

Acknowledgements: I am grateful to O. Ragnisco for pointing out refs. [15,16,17] to me, after 
a first circulation of this paper. 
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